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We report on our recent attempts to determine the short distance behaviors of general 
2-baryon and 3-baryon forces, which are defined from the Nambu-Bethe-Salpeter(NBS) wave 
qq " function, by using the operator product expansion and a renormalization group analysis in 

| QCD. We have found that the repulsion at short distance increases as the number of valence 

quarks increases or when the number of different flavors involved decreases. This global 
tendency suggests a Pauli suppression principle among quark fields at work. 
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§1. Motivation 



The most fundamental quantity in nuclear physics is a force between nucleons, 
the nuclear force, from which major properties of nuclei are extracted. The phe- 
nomenological nuclear potentials P which well describe nucleon-nucleon (NN) scat- 
tering at low energy, exhibit long-to- medium distance attractions, essential for the 
binding of atomic nuclei, as well as short distance repulsion (a repulsive core), im- 
portant for the stability of nuclei against collapse. While the former properties have 
been explained by meson exchanges between nucleons,® the origin of the repulsive 
core® has not yet been well understood. Considerations of the fundamental degrees 
of freedom for nucleons, quarks and gluons, and their dynamics, QCD, are required 
to explain the properties between nucleons at such short distances. 

The recent observation of a neutron star as heavy as twice a solar mass® forces 
us to consider more general short distance repulsions among baryons, not only NN 
^ ■ but also 3N, BB and 3B, where B represents an octet baryon which may include 

one or more strange quarks. This is because although the repulsive core of NN 
potentials is an important ingredient for the determination of the maximum mass of 
neutron stars, it alone seems insufficient to sustain the two-solar-mass neutron star. 
Moreover, an appearance of heavier strange quarks, converted from lighter up and 
down quarks via the weak interaction in such a high density environment, makes the 
equation of state softer in the core of the neutron star, so that the maximum mass 
of the neutron star is further reduced. It has already been pointed out that general 
BB potentials alone can not sustain the two-solar-mass neutron star.^ It is also 
argued that the maximum mass of neutron stars can be increased if 3-body forces, 
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becoming important at high density, also have strong short distance repulsions to 
compensate the reduction due to the appearance of strange quarks. Furthermore it 
is thought that such repulsions should be universal: they should appear not only in 
3NF (3 nucleon forces) but also in 3BF (3 baryon forces)!^' A problem of such an 
explanation, however, is that both experimental and theoretical determinations of 
3-body forces have been known to be very difficult at short distances. 

Recently, a new method has been proposed to define and investigate the NN 
potential from the Nambu-Bethe-Salpeter(NBS) wave function in lattice QCD. This 
method was employed successfully to calculate NN potentials''^ 1 ® and it has been 
widely applied to various hadronic interactions® ffiDEJiGSnSJiEl including BB po- 
tentials at A f v = 3,®® and 3N forces at iV f v = 2,® where A f v denotes the number 
of different quark flavors of the valence quarks of the baryonic systems under con- 
sideration. 

In this paper, we report our recent attempts to determine the short distance 
behaviors of general 2-baryon and 3-baryon forces, using the same definition of po- 
tentials as employed in the lattice QCD simulations referred to above. In our study 
we combine the operator product expansion with a renormalization group (RG) anal- 
ysis of QCD using perturbatively computed RG functions. Our results may give not 
only useful boundary conditions of 2- and 3-body forces at short distance but also 
a hint toward an understanding of origins for short distance repulsions. We have 
observed a tendency toward more short distance repulsions for more nucleons: In 
the case of the NN potential, the short distance repulsion appears only if the ratio 
of two matrix elements has a positive sign.®'^ Although this condition is satisfied 
at low energy,^ the short distance repulsion in not universal. On the other hand, 
a universal short distance repulsion is shown to appear in 3N forces.® We also 
found, for a fixed number of valence quarks, a weaker short distance repulsion for 
increasing N¥ . While the NN potentials (and potentials of general BB systems with 
NJ = 2) have the short distance repulsion at low energy, there appear short distance 
attractions in some channels of BB interactions with = 3PD This tendency is 
found to be true also in 3B forces. 

In the following, we explain how we obtain the above results and what the above 
statements mean more precisely. In sect. El we first give our definition of potentials 
in QCD, which has been employed successfully in lattice QCD calculations. We then 
formulate our analysis methods, involving the operator product expansion and the 
renormalization group in the framework of perturbative QCD. In sect. Owe present 
our results, which include both BB and 3B potentials at both = 2 and 3. A 
conclusion of this paper is given in sect. HI 

§2. Basic idea and methods 

Since potentials themselves are scheme-dependent quantities like running cou- 
plings of QCD, we first give our definition (or scheme) for them. Let us consider the 
equal-time NBS wave function,® which in the case of a BB system for example is 
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defined by 

ip E (f) = (0\B(r, 0)5(0)166, E) in , (2-1) 

where |BB,i±7)i n is a QCD asymptotic in-state with energy E and B{x) is a baryon 
interpolating operator made of 3 quarks such as B(x) = e abc q a (x)q b (x)q c (x) (flavor 
and spin labels suppressed). Although a choice of B{x) is rather arbitrary as long as 
it is coupled to a one baryon state, we take our choice to determine a scheme of the 
NBS wave function (like a choice of the renormalization scheme). It is important to 
note that the N6S wave function contains the information about the BB scattering 
phase shift in QCD at energy Ef^'^ Therefore, from this N6S wave function, the 
non-local but energy-independent 66 potential is naturally defined so as to satisfy 

(e - H )<p E (r) = [ d 3 r'U(ry)cp E (7) , e = — , H = -— , (2-2) 
J m B m B 

for all E < E'th, where E = 2^J k 2 + m B with a baryon mass m# and the threshold 
energy E t ^ = 2ms + tum with a meson mass mu- It is easy to show that the 
non-local but energy-independent U(f, r 1 ) existsP Indeed 

U(?J)= Yl (e-Ho^Eir)^- 1 )^'^) (2-3) 

E,E'<E th 

satisfies eq. (|2-2|) for all E < E t h, where (rj~ 1 ) EE ' is an inverse of 

Vee' = J d 3 np E (r)ip E <(f) = (ifE,<PE>) (2-4) 

in the restricted space where E,E' < Eth- One can construct (r]~ 1 ) EE from eigen- 
vectors corresponding to nonzero eigenvalues of the hermitian operator t]ee'- Note 
that the non-local but energy-independent potential which satisfies eq. (|2-2p for all 
E < £th is not unique, since one can add terms involving wave functions for E > Eth 
which are orthogonal to all wave functions with E < E^. 

In practice, the non-local potential U(r, r 1 ) is expanded in terms of the velocity 
as U(r, r") = V(r, V)5(r — r 1 ), whose lowest few orders for the NN case are given by 

V(r, V) = Voir) + V a {r)a x ■ a 2 + V T (r)S 12 + V LS (r)L • S + 0(V 2 ) , (2-5) 

where r = \r\, &i is the Pauli-matrix acting on the spin index of the i-th nucleon, 
S = (<7i + f?2)/2 is the total spin, L = r x p is the angular momentum, and 
is the tensor operator S\2 = 3(r • <7\)(f ■ a\)/r 2 — o\ ■ This definition of the 
potential has been employed successfully to calculate NN potentials in lattice QCD 
for the first timef^'^'® where the central potentials at the leading order have qual- 
itatively reproduced common features of phenomenological NN potentials: the force 
is attractive at medium to long distance while it has a short distance repulsion, the 
repulsive core. After this success, the method has been widely applied to various 
hadronic mteractionP'^^IH^G! including 66 potentials at N f v = £BM an d 
3N forces at iV f v = 2.Q2J See refPS'ES} for reviews of recent activities. 
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In refp3 we made some toy model studies in order to better understand ques- 
tions related to the scheme dependence, energy dependence, etc. of the NBS wave 
functions and potentials. We noticed that in these 2-dimensional integrable models 
the short distance behavior of the NBS wave function is represented very well by 
the leading terms appearing in the operator product expansion (OPE). The hope is 
that the same qualitative feature applies to QCD, for which, thanks to the property 
of asymptotic freedom, the form of leading short distance behavior of the coefficient 
functions appearing in the OPE can be computed using perturbation theory also for 
nucleons (and baryons). In this section we give a very short summary of our method 
using the example of the NN wave function and potential. Our results for the cases 
NNpSJ BBpD 3]\P3 and 3B will be presented in the next section. 

The behavior of the wave functions i -Pe{'t) at short distances (r = |r| — > 0) is 
encoded in the operator product expansion (OPE) of the two baryon operators: 

B(r/2, 0) B(-r/2, 0) « £ D k (r) O k (0, 0) , (2-6) 

k 

where {Ok} is a set of local color singlet 6-quark operators with two-baryon quantum 
numbers. It is important to note that asymptotically the r-dependence and energy 
dependence of the wave function is factorized into 

<p E (f) «^Z) fc (fO(0|O fe (0,0)|BB,£) in . (2-7) 

k 

Now a standard renormalization group (RG) analysis gives^ 1 the leading short 
distance behavior of the OPE coefficient function as 

D k {r) « A(r)-"* d h , (2-8) 

where A(r) is the 2-loop running coupling defined by 

1 r 

- + ftlnA = In— , r^r*, (2-9) 
A r 

where (in the case of three light dynamical flavors) k = || and r* is some typical 
non-perturbative QCD radius (presumably O(0. 5fm) but a precise specification is 
not needed for our present considerations). 

In (|2-8p v k is related to the 1-loop coefficient of the anomalous dimension of the 
operator O k , d k is the tree- level contribution of 1)^(0). Clearly the operator with 
largest RG power v k dominates the wave function ([2-ip at short distances. We will 
denote the leading (largest) power by v\ and the subleading one (second largest) by 

If neither v\ nor d\ vanishes, this leads to the leading asymptotics of the s-wave 
potential of the form 

27T^T' t 2 ' 10 ) 

m jv r z (In -p-) 

which is attractive for v\ > and repulsive for v\ < 0. Note that although the 
running coupling itself is scheme dependent, the above asymptotic form is not (and 
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it is also energy independent). Scheme dependence only affects the sub-leading 
contributions. Whether (|2- 10[) corresponds to an attractive or repulsive potential is 
easily seen from the formula but can also be understood intuitively from the short 
distance behavior of the corresponding wave function. Namely, if v\ > 0, the wave 
function diverges in the origin (attraction), whereas for v\ < the wave function 
vanishes at the origin (repulsion). 

If v\ = 0, the situation is more complicated. In this case the relative sign of the 
ratio R = (0|O 2 (0, 0)|NN, £) in /(O|0i(0, 0)|NN, E)- m between the leading contribu- 
tion and the subleading contribution corresponding to an operator 2 with v 2 < is 
important. As we will see this occurs in the NN case. If R is positive, the potential is 
repulsive, while it is attractive for negative R. The above intuitive argument seems 
to work also here, namely, for R positive the wave function decreases approaching 
the origin (repulsion) and the other way round for negative R. 

To determine the asymptotic behavior we thus have to compute the spectrum 
of the renormalization group power corresponding to the operators appearing in 
the BB OPE. The general form of a gauge invariant local 3-quark operator is given 
by 

B F r {x) = B f J*(x) = e ahc q a J{x) q f{x)q^{x) , (2-11) 

where a, j3, 7 are spinor, /, g, h are flavor, a, b, c are color indices of the (renormalized) 
quark field q. The color index runs from 1 to N c = 3, the spinor index from 1 to 4, 
and the flavor index from 1 to . As indicated in (|2-lip we use capital labels for 
sets of flavors F = fgh and Dirac labels r = a/3j. Note that B^ is symmetric 

under any interchange of pairs of indices (e.g. B*JL = B 9 ^ ) because the quark fields 
anticommute. The usual nucleon operator which is employed in lattice simulations 
is a certain linear combination constructed from the above operators.^ 

Gauge invariant local 6-quark operators can be written as linear combinations of 
simple operators = Bp* B r \ These local operators mix among themselves 

under renormalization. The one-loop mixing matrix is given by diagrams corre- 
sponding to the exchange of one gluon between all pairs of quark lines. Due to the 
flavor symmetry (and partially to the chiral symmetry) of massless QCD, only those 
6-quark operators (of the form [BB\p^p^) which preserve the set of flavor and Dirac 
indices as 

f 1 uf 2 = f a uf b , r 1 ur 2 = r A ur B , (2-12) 

occur in the mixing problem. Note however that such operators are not all linearly 
independent. We have the following "gauge" constraints F^l 

3[^SS+E[5B]g^=0. (2-13) 

Here means a simultaneous exchange between the i-th indices of Fi,Ti and the 
j-th indices of ^2,^2- 

*' We call these constraints "gauge" constraints because we used them to show the gauge in- 
variance of the 1-loop anomalous dimensions. 
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The problem is thus reduced to finding the set of independent operators af- 
ter imposing the gauge constraints and then writing down and diagonalizing the 
corresponding one-loop mixing matrix. This task can easily be automated and we 
calculated and checked the results using independent Mathematica and Maple pro- 
grams. The size of the linear algebra problem turned out quite large in some cases. 
For example the largest dimensions occur for a certain set of Dirac indices in the 
BBB case: there are originally 14130 operators before and 1518 after imposing the 
constraints. We are indebted to Thomas Hahn for providing us with effective linear 
algebra tools for solving this part of the problem with Mathematica. 



§3. Results 



In the first three subsections we briefly summarize our results for the leading 
behavior of 2 nucleon, 2 baryon (involving operators with 3 different flavors) and 

3 nucleons; for the details of these computations the reader is requested to con- 
sult refs.^'^'^ respectively. Our latest results concerning the case of 3 baryon 
operators are presented in the fourth subsection. 

In each case the computation consists of three steps, i) a basis of independent 
gauge invariant operators of engineering dimension equal to that of the operator 
product appearing in the definition of the BS wave function is determined for each 
Dirac structure, ii) The spectrum of 1-loop anomalous dimensions 7& of the linear 
combination of these operators which do not mix under renormalization (at 1-loop) 
are computed, iii) The operators appearing in the OPE at tree level are listed. 

For operators which renormalize multiplicatively (without mixing) 

O ren = Z (e,<7)O w , (3-1) 

the associated anomalous dimensions are, using dimensional regularization in D = 

4 — 2e dimensions, specified by 

TO = - [eg + /3 05 3 + • • • ] ^ = 7cV + 0(<? 4 ) , (3-2) 



where /3q is the QCD 1-loop beta function coefficient 

(3-3) 



ft ' 



16vr 2 



11 - -N { 

3 1 



We will restrict attention to the octet baryon operators appearing in the definition 
of the BS wave functions, which have 1-loop anomalous dimension 

7B0 = 24d, d = l/(96vr 2 ). (3-4) 

It turns out that all 1-loop anomalous dimensions of operators appearing below can 
be expressed as even integers times 2d. 

3.1. The NN potentials 

The computationally easiest case is the BS wave function defined with the prod- 
uct of two nucleons, operators defined with just two (u, d) flavors of quarks. The 
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spectrum of anomalous dimensions in this case has a striking structure: 

7fco<27 jBO , Vfc. (3-5) 

It is easy to see that there are NN (also BB) operators with 7^0 = 27^0 - e.g. those 
of the form B + B^ where B± involves only quarks of positive (negative) chirality. 
Such operators generically appear in the OPE at tree level and hence in the NN 
case these dominate the short distance behavior of the BS wave function. However 
since = const, asymptotically for these operators, their contributions alone do 
not determine the short distance behavior of the potentials (defined through taking 
derivatives of the wave function). The leading behavior of the (central) potentials 
involves also the next-to-largest anomalous dimensions, which are found to be 7( 0,1 ) = 
24d , 7 (1,0) = 40d in the two spin-isospin channels (S,I) = (0, 1), (1,0) respectively: 

K c (S ' J) (r) ~ R^r-^-lnr/ 8 ^- 1 , (3-6) 
with (assuming Nf is such that /3o > 0) 

^W- 3^LzjM <0 , (3-7) 
2po 

and the constants R^ s '^ appearing above involve the ratio of matrix elements of the 
leading and subleading operators in the OPE as 

hbo = (gggWk (3 . 8) 

(0|oP' , |NN,E) il , 

Thus although the OPE determines the functional form of the short distance behav- 
ior, unfortunately it does not determine the sign in this special case. Information 
on the sign of R&D requires additional non-perturbative input, which could e.g. be 
supplied by dedicated numerical simulations. In^ ;) we have given arguments based on 
a non-relativistic expansion that the are positive, in which case the potential 

is repulsive in both channels as indeed observed in the numerical simulations. 

3.2. The BB potentials 

We next consider the operator product of two baryon octet operators where the 
operators appearing in the OPE can involve 3 different flavors. The tensor product 
of two octets can be decomposed into 6 channels, according to 

8 8)8= (10 8 27), + (80 10 ©IU) a , (3-9) 

where subscripts s(a) represents the property that the operator is (anti-)symmetric 
under the interchange of two baryons. The two nucleon states with (S, I) = (0, 1) or 
(1,0) considered above belong to the 27 s or 10 a representations, respectively. It is 
sufficient to consider the situation of composite 6-quark operators with two quarks of 
each flavor (uuddss) since all 6 SU(3) representations appear in this type of operator. 

Although most of the operators have 1-loop anomalous dimensions < 27^0 we 
found also operators with anomalous dimensions greater than 27^0 corresponding 



8 



S. Aoki, J. Balog, P. Weisz, 



Table I. List of channels with anomalous dimensions greater than 2-^bo\ others can be obtained by 
symmetry transforms of the Dirac indices lf>2,3H4or (1)2) o (3, 4). 



Dirac structure 


7S7(2d) 


SU(3) representation 


112334 


42 


1 


111222 


36 


1 


112234 


36 


1 


111234 


32 


108 


112234 


32 


108 



to attractive BS potentials in the three channels l s ,8 a ,8 s ; they are listed in Table HI 
There are no attractive operators in the 27 s ,10 a ,10 a channels, which is consistent 
with what we found in the two nucleon case. All the attractive operators appear in 
the OPE of two baryon operators at tree level and therefore the corresponding BS 
potentials in these channels (denoted by 1Z) behave asymptotically as in eq. (|2-1U|) 
with 

>0 . (340) 

The larger anomalous dimensions occur in the lower dimensional SU(3) representa- 
tions. Indeed the tendency of repulsion among the operators with fixed total number 
of valence quarks is increased if the number of participating valence flavors is de- 
creased; this is a first indication of what we suggestively call a "Pauli suppression 
principle" at work. 

Monte Carlo simulations in 3-flavor lattice QCEP^'^ 1 have as yet only found 
attraction in the l s channel. The fact that attraction in octet channels has not yet 
been observed may have many different reasons. Firstly we note that the amplitudes 
of attractive operators in the 8 S channel vanish in the non-relativistic limit. Secondly, 
and this comment applies to all our considerations, we do not know at which distance 
r the asymptotic formulae (|2-10p set in. The physical distances of the BS potentials 
probed in Monte Carlo simulations up to now (>~ 0.2 fm) may not be small enough. 

3.3. The 3N potentials 

For the investigation of the 3 baryon potentials we consider equal-time NBS 
wave functions involving the product of three baryon operators: 

V>3B,Mr,p) = (0\B(£i, 0)B(x 2 , 0)B(x 3 , 0)|3B,£) in , (3-11) 

where \SB,E) i n is a 3-baryon asymptotic in-state of energy E. In the argument of 
ip3B,E we have introduced Jacobi coordinates r = x\—xi ,p= [x 3 — (%i + ^2)/2] /\/3. 
The three baryon potential is defined through the wave functions by 



— (v r 2 + vj) + V 3B 

TUB 



V 3B (r,p)= VBBtfi-fy + VsBftp), (3-13) 

1<«<7<3 
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where the reduced masses of these Jacobi coordinates are m^/2 and only the leading 
contributions in the velocity expansion are presented. The Vbb are the BB potentials 
defined through the BB wave functions. 

The leading terms in the OPE of 3-baryon operators in r , p — > limits involves 
the associated local gauge invariant 9-quark operators (of minimal engineering di- 
mension of 27/2). 

For the 3N case the spectrum of anomalous dimensions of the the 9 quark oper- 
ators (now involving only two flavors) have the simple property: 

7ao<3 7 bo, VA. (3-14) 

Again it turns out that most of the operators with largest anomalous dimensions 
appear in the OPE at tree level. From this result, without recourse to additional 
nonperturbative arguments, we deduce in this case that generically the total central 
3N potential V^{r,p) is repulsive at short distances: 

V 3 N(r,p)~ 2 4 f^ (3-15) 



where s = y (r) 2 + (p) 2 and 

Pa = 7A ° no 7B ° < °) max7Ao = 16(i. (3-16) 

Since the singularity of Vnn given in (|3-6|) is milder than that of V3N in (|3T5|) . it 
follows that the 3N potential V3N defined in (|3T3p (with B — > N) is responsible for 
the asymptotic repulsion. 

We note that increasing the number of quarks, but not the number of flavors, 
tends to decrease the anomalous dimensions of the operators, thereby increasing the 
tendency for a repulsive core of the potential, which is a further indication of a "Pauli 
suppression principle" mentioned before. 

3.4. The 3B potentials 

This problem is very similar to the previous three cases but since the corre- 
sponding results have not yet been published we now describe the outcome of our 
considerations in some detail. However, due to the very large dimension of the prob- 
lem we refrain from giving full details in the paper. The tensor product of three 
octets can be decomposed as 

8 © 8 <g> 8 = 64 © (35 © 35) 2 © 27 6 © (10 © TU) 4 © 8 8 © 1 2 (3-17) 

where the subscript indicates the multiplicity of the corresponding irreducible rep- 
resentation. 

In the three-baryon case we are considering 9-quark operators with three quarks 
of each flavor (uuudddsss). We find that there are 35 different Dirac structures (plus 
their transforms under the symmetries 1 -H- 2, 3 -H- 4 or (1, 2) -H- (3, 4)). The simplest 
Dirac structure is 111111111 and there are 9 such operators but obviously only one of 
them (corresponding to ria=i(' u i^i s i)) remains after imposing the gauge constraints. 
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Table II. The Dirac case 111111122. 



7V(2d) 


SU(3) representation 


-52 


27 


-36 


1,8 


-22 


8 


-4 


1 



The anomalous dimension of this operator is —36 (in 2d units) and is an SU(3) singlet. 
The next simplest case corresponds to the Dirac structure 111111112 where there 
are originally 51 operators but only 3 remain independent after imposing the gauge 
constraints. The remaining three operators correspond to anomalous dimension —36 
and form a singlet and the 1 = and 1 = 1 members of an SU(3) octet. Next 
we consider the Dirac case 111111122. Here the number of operators is reduced 
from 153 to 9. The results are summarized in Table HI] For more complicated Dirac 
structures the number of operators increases rapidly and for example for 111333224 
and 111222334 we have 9960 operators, which are reduced to 1014 by the gauge 
constraints. The largest case is 111223344. The number of operators here is reduced 
from 14130 to 1518. The spectrum of these 1518 independent operators is 
{-76, -60, -54, -52, -48, -46, -44, -42, -40, -36, -34, 
-30, -28, -24, -22, -20, -18, -16, -12, -10, -8, -6, -4, 
0, 2, 4, 6, 8, 12, 14, 18, 20, 24, 26, 30, 32, 36, 38, 42, 44, 48} 
and the dimensions of the SU(3) representations in which these operators transform 
are 

{1, 8, 10, 27, 28, 35, 55, 64, 80, 81} . 
All eigenvalues and the corresponding SU(3) representations are given in Table IIVI 
for this case. 

The most important new feature of the spectrum of operators with three-baryon 
quantum numbers (with respect to the three-nucleon case) is the presence of at- 
tractive channels, though the vast majority of eigenvalues correspond to repulsion 
as before. We give a complete list of such channels (corresponding to eigenvalues 
7A0 > 37so) hi Table ITTTl where an asterisk superscript (*) indicates if the operator 
is already present at tree level in the product of three baryon operators. 

As can be seen from the table, we have found tree operators with 7^0 > hi 
the 1 (singlet) and 8 (octet) channels, which give a universal attraction at short dis- 
tances (attractive cores). In addition, there are tree level operators with 7^0 = 37bo 
in the 8 (octet) and 10, 10 (decuplet) channels, which lead to non-universal attrac- 
tion or repulsion at short distances depending on the sign of the ratio of associated 
matrix elements. On the other hand, universal repulsion at short distances (repulsive 
core) appears only in the 64, 35, 35, 27 channels. 

Note however, that the anomalous dimensions listed in the table give, using 
(|3-15|) and (|3T6[) . the asymptotic behavior of the "total" 3B potential V3B and not 
the "true" 3B potential V3B, which is obtained from the former by subtracting the 
sum of the 2B potentials. In the 3N case this distinction was not important, since the 
leading behavior of the 2N potentials is always milder than that of the 3N potential, 
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Table III. List of channels with anomalous dimensions greater than or equal to 3^bo- Operators 
present at tree level in the OPE are indicated by an asterisk *. 



±y 1.1 oil Llt_< Li hi c 




T ( *^ 1 TTTIT'PQPTlt" i1"t"H~lt"l 




48 


1 8* 
1, 8 




AA 


1*8* 




AO 


1*8* 




oo 


1 8 
1, 8 




36 


8*, 10*, TO* 


111333224 


44 


r,s* 




38 


1,8 




36 


1 


111222334 


48 


8* 




44 


1* 


111133442 


42 


1*,8* 




38 


1,8 


111122334 


44 


1* 




36 


1 


111133324 


38 


1,8 


111113324 


36 


1 



which is hence the dominant contribution at short distances. The 3B case is made 
more complicated by the fact that both the "total" 3B potential V3B and the 2B 
potential V2B have some attractive channels, competing with each other. 

Take for example the 3B configuration with Dirac structure 111223344, which 
contains the largest anomalous dimension 48 (in 2d units). This is in the octet 
channel and leads to an attractive total force, whose strength is proportional to the 
numerator of the formula (|316[) . 12 in our units. On the other hand, there is a 
2B contribution to this total force coming from two baryons with Dirac structure 
112334 and forming a singlet state with the highest anomalous dimension 42 in our 
units (times a free octet baryon). From (|3T0p we see that the corresponding 2B 
force is proportional to 18. So in this example the 2-body attractive potential is 
actually stronger than the total 3-body attraction hence there must be a repulsive 
effect coming from V3B as 



V 3 B(f,p) ~ — 

m B 



s 



- ln(s/r,)) rjJ-(- H r ij/ r *)) 



(3-18) 



as \r\, \p\ — > 0, where mg is the octet baryon masa^], and ( for Nj = 3) 



2 „ 1 



/3 3 m B X =g, /3g 1 f| x = 3, s = V^+W, (3-19) 
T12 = M, r 13 = \f/2 - V3p\, r 23 = \f/2 + V3p\. (3-20) 

As an other similar example we take the Dirac structure 111133442 in the singlet 
channel with anomalous dimension 42. A two-baryon operator with Dirac structure 



We here ignore mass differences among octet baryons caused by the flavor SU(3) breaking. 
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111234 in the octet channel (combined with an extra octet baryon) contributes to 
this component of the 3B force and the corresponding numbers are: 

/3 3 m B ax = ^ /3f?B ax = ^- (3-21) 

We do not further attempt to analyze the general case, which is clearly very 
complex since in addition to the cancellation between a 2B and a 3B force as in the 
above examples, there are cases where cancellation occurs among the 2B forces. 

§4. Conclusion 

In this paper, we have presented our recent activities on determinations of short 
distance behaviors of BB and 3B potentials defined from the NBS wave function 
in QCD using the operator product expansion and the RG analysis in perturbative 
QCD. Our results show that the repulsion at short distance (repulsive core) become 
stronger for more participating valence quarks or less different numbers of flavors 
NJ . Explicitly we have: 

1. NN potentials (iV7 = 2) seem to have a repulsive core at low energy, which 
however is not universal in the sense that the coefficient may depend on the 
properties of the NBS wave function such as total energy. 

2. BB potentials (N^ = 3) can have not only a repulsive core but also channels with 
an attractive core, the latter of which is universal (i.e. energy independent). 
In particular, the attractive core in the flavor singlet potential has indeed been 
confirmed in lattice QCD simulations.^'^ 

3. 3N potentials have a universal repulsive core. This is the most unambiguous 
result in our project. 

4. In the 3B case we were able to study the short distance asymptotics of the 
"total" 3B potential V3B only with our methods. Although the vast majority 
of various channels (64,35,35,27) in 3B potentials have a universal repulsive 
core, a few channels (1,8) have a universal attractive core. In addition, there 
may appear a non-universal repulsive or attractive core in the 10, 10 channels. 

Roughly speaking, these results suggest that a "Pauli suppression principle" among 
quarks is at work. Explicit 1-loop calculations are needed, however, to obtain the 
detailed structure of the anomalous dimensions of BB or 3B operators. 

Our results give important information on the short distance behaviors of bary- 
onic interactions, which can be used to constrain baryonic potentials at short distance 
obtained in lattice QCDp9»23 in particular, to constrain 3B potentials-^- . As men- 
tioned in the introduction, the universal repulsion of 3B potentials seems necessary 
to explain an existence of two-solar-mass neutron star. Our result shows that the 
vast majority of 3B channels have a repulsive core in the 3B potential, though a few 
channels have an attractive core. It is interesting and important to check these qual- 
itative behaviors of 3B potentials at short distances and furthermore to investigate 
the strength of these cores, employing lattice QCD simulations. 
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Table IV. The Dirac case 111223344. Operators present at tree level in the OPE are indicated by 
an asterisk *. 



7o7(2d) 


SU(3) representation (multiplicity) 


48 
44 


1,8* 
1*,8* 


42 


1*,8* 


38 


1,8 


36 


8*, 10* 


32 


1*,8*(3),10*(2) 


30 


8, 10* 


26 


1,8(3), 10(2) 


24 


1,8*(4),10*(3),27*(2) 


20 


1*(2),8*(6),10*(4),27*(3) 


18 


1,8*(4),10*(3),27*(2) 


14 


1(2), 8(5), 10(3), 27(3) 


12 


1(3),8*(7),10*(5),27*(4),35* 


8 


1,8*(8),10*(9),27*(5),35*(2) 


6 


1*(2),8*(6),10*(6),27*(4),35*(2) 


4 


1,8(3), 10(2), 27 


2 


8(6), 10(9), 27(5), 35(3) 





1(3), 8*(9), 10* (10), 27* (9), 35* (5), 64 


-4 


1* (4), 8* (9), 10* (10), 27* (9), 35* (6), 64 


-6 


1(2), 8* (8), 10* (9), 27*(7), 35(4), 64 


-8 


8(2), 10(3), 27, 35 


-10 


1(2), 8(7), 10(9), 27(8), 35(6), 64(2) 


-12 


1,8(5), 10(8), 27(5), 35(5), 64 


-16 


1(3), 8* (9), 10* (11), 27* (11), 28, 35* (9), 64(3) 


-18 


1, 8* (7), 10* (12) , 27* (9), 28, 35* (10), 64(3) , 81 


-20 


1,8(4), 10(5), 27(4), 35(3), 64 


-22 


1, 8(5), 10(9), 27(8), 28(2), 35(9), 64(2), 81 


-24 


1, 8* (4), 10(5), 27* (7), 28(4), 35* (9), 64(3), 81 


-28 


1(2), 8(6), 10(7), 27(7), 28(2), 35(9), 64(4), 81(2) 


-30 


1, 8(4), 10(6), 27(7), 28(2), 35(8), 64(3), 81 


-34 


1, 8(5), 10(6), 27(6), 28, 35(7), 64(4), 81(2) 


-36 


8(3), 10(4), 27(4), 28, 35(5), 64(3), 81(2) 


-40 


8, 10(3), 27(3), 28(2), 35(5), 64, 81(2) 


-42 


8(2), 10(3), 27(3), 28, 35(4), 64(2), 81(2) 


-44 


8, 10,27,35,64,81 


-46 


8, 10(3), 27(3), 28(2), 35(5), 64(2), 80, 81(2) 


-48 


27, 28,35,64,80,81 


-52 


10, 27(2), 28(2), 35(3), 64(2), 80(2), 81(2) 


-54 


10, 27, 28, 35(2), 64, 80, 81 


-60 


10, 27(2), 28(2), 35(3), 64(2), 80(2), 81(2) 


-76 


35,55,64,80,81 



